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ABSTRACT. We take some first steps in providing a synthetic theory of dis- 
tributions. In particular, we are interested in the use of distribution theory as 
foundation, not just as tool, in the study of the wave equation. 
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Introduction 

The aim of this paper is to contribute to a synthetic theory of distributions. 
The sense in which we understand "synthetic" in this context is that we 
place ourselves in a setting (category) where everything is smooth (differ- 
entiable). Now distributions are sometimes thought of as very non-smooth 
functions, like the Heaviside function, or the Dirac delta. We take the view- 
point, stressed by Lawvere, that distributions are extensive quantities, where 
functions are intensive ones. It is only by a spurious comparison with func- 
tions that distributions seem non-smooth. 

A main assumption about the category in which we work is that it is 
cartesian closed, meaning that function- "spaces", and hence some of the 
methods of functional analysis, are available. 

This viewpoint also makes it quite natural to formulate the wave equation 
as an evolution equation, i.e. an ordinary differential equation describing the 
evolution over time of any initial distribution, so it is an ordinary differential 
equation with values in the vector space of distributions. 

The main construction in the elementary theory of the wave equation is 
the construction of the fundamental solution: the description of the evolution 
of a point (Dirac-) distribution over time. (Other solutions with other initial 
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states may then by obtained by convolution of the given initial state with 
the fundamental solution; we shall not go here into this classical technique.) 

To say that distributions are extensive quantities implies that they trans- 
form covariantly. To say that functions are intensive quantities implies that 
they transform contravariantly. Distributions are here construed, following 
Schwartz, as linear functionals on the space of (smooth) functions. But since 
all functions in the synthetic context are smooth, as well as continuous, there 
is no distinction between distributions and Radon measures. 

So we consider a cartesian closed category E with finite limits, in which 
there is given a commutative ring object R, to be thought of as the real 
number line. 

Already on this basis, one can define the vector space T>' C (M) of distribu- 
tions of compact support on M, for each object M G E, namely the object 
of i?-linear maps R — > R ("vector space" in this context means i?-module). 

We shall assume that elementary differential calculus for functions R — > R 
is available, as in all models of SDG, cf. [I], [TT], [Hj, etc. We shall also assume 
some integral calculus, but only in the weakest possible sense, namely we 
assume 

Integration Axiom: For every ip : R — > R, there is a unique \I> : R —>■ R 
with \P' = if) and with *(0) =0. 

Note that we do not assume any order < on R, so that "intervals" 
[a, b] C R do not make sense as subsets. "Intervals" , on the contrary, will be 
construed as distributions: for a,b G R, [a, b] denotes the distribution 



The right hand side here of course means ty(b)—^(a), where *Sf is the primitive 
of if; given by the integration axiom. (This weak form of integration axiom 
holds in some of the very simple models of SDG, like in the topos classifying 
the theory of commutative rings.) 

Finally, for the specific treatment of the wave equation, we need that 
the trigonometric functions cos and sin should be present. We assume that 
they are given as part of the data, and that they satisfy cos 2 + sin 2 = 1, 
and cos' = — sin, sin' = cos. Also as part of the data, we need specified an 
element it G R so that cos7r = — 1, cosO = 1. 
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Except for the Taylor Series/ Formal Solutions considerations in the end 
of the paper, the setting does not depend on the "nilpotent infinitesimals" of 
SDG, but could also be, say, that of Froelicher-Kriegl [2 , or Grothendieck's 
"Smooth Topos". 

We would also like to remark that one can probably construct such smooth 
toposes in which no non-trivial distribution of compact support has a den- 
sity function, or equivalently, no function (other than 0) gives rise to a 
distribution-of-compact-support; our description of fundamental solutions to 
the wave equation would not be affected. An example seems to be the topos 
classifying C^-algebras, where C u is the algebraic theory of entire real- or 
complex- analytic functions. 

1 Generalities on distributions 

We want to apply parts of the general theory of ordinary differential equations 
to some of the basic equations of mathematical physics, the wave- and heat- 
equations 1 . 

This takes us by necessity to the realm of distributions. Not primarily 
as a technique, but because of the nature of these equations: they model 
evolution through time of (say) a heat distribution. A distribution is an 
extensive quantity, and does not necessarily have a density function, which 
is an intensive quantity; the most important of all distributions, the point 
distributions (or Dirac distributions), for instance, do not. 

As stressed by Lawvere in [SJ, functions are intensive quantities, and 
transform contravariantly; distributions are extensive quantities and trans- 
form covariantly. For functions, this is the fact that the "space" of functions 
on M, R M is contravariant in M, by elementary cartesian-closed category 
theory. Similarly, the "space" of distributions of compact support on M is a 
subspace of R R (carved out by the R linearity condition), and so for similar 
elementary reasons is covariant in M. 

Let us make the formula for covariant functorality T>' c explicit. Let / : 
M -> iV be a map. The map V' c (f) : V' C (M) -> V' C {N) - which may also be 
denoted /* - is described by declaring 

< f*(T ), <P >=< T, o / >, (1) 

1 We do not discuss the heat equation in the present paper; we hope to return to it and 
improve the version of [Hj 
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where T is a distribution on M, and is a function on N. The brackets 
denote evaluation of distributions on functions. If we similarly denote the 
value of the contravariant functor M i— > R M on a map / by /*, the defining 
equation for /* goes < f*(T), >=< T, f*(<f>) >■ 

We note that T>' C (M) is an i?-linear space, and all maps /* : V' C (M) — > 
T>' C {N) are i?-linear. Also V' C (M) is a Euclidean vector space V, meaning that 
the basic differential calculus in available, for instance that the basic axiom 
of SDG holds; we return to this in Section 2. 

For any distribution T of compact support on M, one has its Total, which 
is just the number < T, 1 >G R, where 1 denotes the function on M with 
constant value 1. Since /*(1) = 1 for any map /, it follows that /* preserves 
Totals. (Alternatively, let 1 denotes the terminal object (=one-point set). 
Since V c (l) = R canonically, the Total of T may also be described as !*(T), 
where ! : M — * 1 is the unique such map. Then preservation of Totals follows 
from functorality and from uniqueness of maps into 1.) 

Recall that a distribution T on M may be multiplied by any function 
g : M — > R, by the recipe 



A basic result in one- variable calculus is "integration by substitution". 
We present it here in pure "distribution" form; note that no assumption on 
monotonicity or even bijectivity of the "substitution" g is made. 

Proposition 1 Given any function g : R — > R, and given a,b € R. Then, 
as distributions on R, 



Proof. Let if) be a test function, and let \l/ be a primitive of it, ^' = if>. 
So < [g(a), g(b)],ip >= ^(g(b)) - V(g(a)). On the other hand, by the chain 
rule, ^ o g is a primitive of g' • (if) o g), and so 



■V(g(a))--*(g(b)) =< [a, 6], g'-tyog) >=< g'-[a,b],ifjog >=< g^g'-^b}),^ > . 



<g-T,<f>>=<T,g-<f>> . 



(2) 



g*{g' ■ [a,b]) 



[g(a),g(b)}. 
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The external product of distributions of compact support is defined as follows. 
If P is a distribution on M, and Q a distribution on N, we get a distribution 

P x Q on M x N, by 

< P x Q,ijj >=< P, [m (->■< Q, -)>]>■ 

In general, the external product construction x will not be the same as the 
external product construction ~x given by 2 

< PxQ,ip >=< Q, [n i->< P,ip(-,n) >] > . 

However, if [a, b] and [c, d] are intervals (viewed as distributions on R, as 
described above), [a,b] x [c,d] = [a,b\x[c,d], as distributions on R 2 , by an 
application of Fubini's Theorem, (which holds in the context here - it is a 
consequence of equality of mixed partial dervatives). - Distributions arising 
in this way on R 2 , we call rectangles. The evident generalization to higher 
dimensions, we call boxes. We have 

pb pd 

< [a, b] x [c, d] , ip >— / / ip(x,y) dy dx, 

J a Jc 

in traditional notation. Notice that we can define the boundary of the box 
[a, b] x [c, d] as the obvious distribution on R 2 , 

(P 2 c )*[a, b] + (pl)*[c, d] - (p 2 d )*[a, b] - (pl)*[c, d] 

where p 2 (x) = (x, c), p\(y) = (b,y), etc. 

By a singular box in an object M, we understand the data of a map 
7 : R 2 — > M and a box [a, b] x [c, d] in R 2 , and similarly for singular intervals 
and singular rectangles. Such a singular box gives rise to a distribution on 
M, namely g*([a,b] x [c,d]). 

By "differential operator" on an object M, we here understand just an R- 
linear map D : R — > R M . If D is such an operator, and T is a distribution 
on M, we define D(T) by 

<D(T ),jj >:=< T,D(ij) >, 

In fact, the two external product formations described here provide the covariant 
functor T>' c (— ) with two structures of monoidal functor E — > E, in fact, they are the 
monoidal structures that arise because T>' c {—) is a strong functor on E with a monad 
structure, 0, 0- 
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and in this way, D becomes a linear operator V C (M) — ► T>' C (M). 

In particular, if X is a vector field on M, one defines the directional 
derivative Dx{T) of a distribution T on M by the formula 

<D x (T),i/>>=<T,D x (if>)> . (3) 

This in particular applies to the vector field d/dx on R, and reads here 
< T',i/j >=< T,ij}' > {ip' denoting the ordinary derivative of the function 
ij}). (This is at odds with the minus sign which is usually put in into the 
definition of T', but it will cause no confusion - we are anyway considering 
second order operators, where there is no discrepancy.) 

The following Proposition is an application of the covariant functorality 
of the functor T>' c , which will be used in connection with the wave equation 
in dimension 2. We consider the (orthogonal) projection p : R 3 —> R 2 onto 
the xy-plane; A denotes the Laplace operator in the relevant R n , so for R 3 , 

A is d 2 /dx 2 + d 2 /dy 2 + d 2 /dz 2 . 

Proposition 2 For any distribution S (of compact support) on R 3 , 

p.(A(S)) = A(p„(fif)). 

(The same result holds for any orthogonal projection p of R n onto any linear 
subspace; the proof is virtually the same, if one uses invariance of A under 
orthogonal transformations.) 

Proof. For any tp : R 2 — > R, 

A(pV))=p*(A(^)), 

namely d 2 ip/dx 2 + d 2 ip/dy 2 . From this, the Proposition follows purely for- 
mally. 

2 Calculus in Euclidean vector spaces 

Recall that a vector space in the present context just means an i?-module. 
A vector space E is called Euclidean if differential and integral calculus for 
functions R — > E is available. An axiomatic account is given in jl], [TI], jS] 
and other places. The coordinate vector spaces are Euclidean, but so are also 
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the vector spaces R M , and V C (M) for any M. To describe for instance the 
("time-") derivative / of a function / : it! — > V C (M), we put 

< >=| </(*), V>>- 

Similarly, from the integration axiom for R, one immediately proves that 
T>' C (M) satisfies the integration axiom, in the sense that for any h : R — > 
V' C (M), there exists a unique H : i? — > T>' C (M) satisfying ii(0) = and 
H'(t) = h(t) for all t. In particular, if /i : R -»• 2^(M), the "integral" 
Jah(u) du makes sense (as H(b) — H(a)), and the Fundamental Theorem of 
Calculus holds, almost by definition. 

As a particular case of special importance, we consider a linear vector 
field on a Euclidean i?-module V. To say that the vector field is linear is to 
say that its principal-part formation V — > V is a linear map, T, say. We have 
then the following version of a classical result. By a formal solution for an 
ordinary differential equation, we mean a solution defined on the set of 
nilpotent elements in R (these form a subgroup of (R, +)). 

Proposition 3 Let a linear vector field on a Euclidean vector space V be 
given by the linear map T : V — > V. Then the unique formal solution of the 
corresponding differential equation, i.e., the equation F(t) = F(F(t)) with 
initial position v, is the map x V — > V given by 

(i, u )_ e *-r( u ) > (4 ) 

where the right hand side here means the sum of the following "series" (which 
has only finitely many non-vanishing terms, since t is assumed nilpotent): 

v + tT(v) + t ^(v) + t ^(v) + --- 
(Here of course T 2 (v) means T(r(f)), etc.) 

Proof. We have to prove that F(t) = T(F(t)). We calculate the left hand 
side by differentiating the series term by term (there are only finitely many 
non-zero terms): 

m + 1 • r» + |^r» + ... = r> + 1 ■ r» + 1 • r» + ...) 
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using linearity of T. But this is just T applied to F(t). 

There is an analogous result for second order differential equations of the 
form F {t) = T(F(t)) (with T linear); the proof is similar and we omit it: 



Proposition 4 The formal solution of this second order differential equation 
F— FF, with initial position v and initial velocity w, is given by 

F( t ) = v + t-w + |jr» + + ^T 2 (v) + |r» + .... 



We shall need the following result ( "change-of- variable Lemma"); for V = R, 
it is identical to Proposition Q and the proof is in any case the same. 

Proposition 5 Given f : R —>V, where V is a Euclidean vector space, and 
given g : R — > R. Then for any a, b £ R, 

/ f(g{x)) ■ g'( x ) dx= f(u) du. 



Linear maps between Euclidean vector spaces preserve differentiation and 
integration of functions R — > V; we shall explicitly need the following par- 
ticular assertion 

Proposition 6 Let F : V — >• W be a linear map between Euclidean vector 
spaces. Then for any f : R — > V , 

F([ b fit) dt)= f b F(f(t))dt 

J a J a 



3 Spheres and balls as distributions 

Let S be a distribution in R n ; ultimately, it will be the unit sphere, see below. 
We describe some families of distributions derived from it. Let t £ R (not 
necessarily t > - we haven't even assumed an order relation on R). We 
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then have the homothety "multiplying by t from R n to R nn , which we denote 
H\ so 

H\x) = t-x, 

for any x G R n . 

We are going to use the covariant functorality of T>' c with respect to these 
maps H l . Note that for any distribution T on R n , 

H°(T) = Total(T)-5(0), (5) 

where <5(0) denotes the Dirac distribution at G R n , given by < 5(0), ip >= 
ip(0) . We put 

S* ■= Ht(S). 

It has the same Total as S, but its support 3 is larger (e.g. for t — 2, it is 
2™" 1 times as big as that of S). So if S is the unit sphere, S l is "the diluted 
sphere of radius t" . We also want an undiluted sphere of radius t; we put 

S t := f 1 - 1 ■ SK 

Note that in dimension 1, S t = S*. 

The ball of radius 1 is made up from shells (undiluted spheres) "of radius 
u (0 < u < 1)" (heuristically !), motivating us to put 

B :— l S u du, 
Jo 

using integration in V' c (R n ). Let t G R. We put 

B l := Hl(B). 

It has the same Total as B, but its support is larger ( "if t > 1" - heuristically), 
so if S is the unit sphere, B is u the diluted ball of radius t" (think of the 
expanding universe). We also want an undiluted ball of radius t; we put 

B t := t n ■ B l . 

We then have 

3 We haven't here introduced the notion of support of a distribution, and only use the 
word here for motivating the word "diluted" . 
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Proposition 7 For all t G R, 

B t = [ S v dv. 
Jo 



Proof. 

B t = t n ■ B f = t n ■ Hl(B) 

= t n -Ht([ 1 S u du) 
Jo 

= t n - f 1 Hl(S u )du 
Jo 

(by Proposition IHJ) 

= t n f 1 Hi{u n - x ■ h:(s)) du 

Jo 

= t n ■ [ l u^H^iS) du 
Jo 

= f 1 (t-u) n - 1 ■ Hl u (S)tdu 
Jo 

= f v n - x Hl{S) dv 
Jo 

(by change-of- variable Lemma (Proposition EJ), with v := t ■ u), which is 
Jq S v dv, as claimed. 

We now give explicit defining formulae for S in dimensions 1, 2 and 3. 
These are of course standard integral formulae in disguise - explicit integral 
formulae come by applying the definitions, and then integral formulae for S l , 
St and B t may be derived (using Proposition [7| and related arguments) - we 
give some of these formulae below. 

Dimension 1 S := 5(1) + S(-l) 

Dimension 2 S := cis*([0, 2tt]), where cis : R — » R 2 is the map 9 h- > 
(cos 9, sin 9). 

Dimension 3 S := s ■ sph J|J ([0, 27r] x [0, 7r]), where s : R 2 — > R is the function 
(9, (ft) i—*- sin 0, and where sph is "the spherical coordinates map" R 2 — ► R 3 
given by 

(9, <j)) i — ► (cos 6* sin 0, sin 6* sin 0, cos 0) . (6) 
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In dimension 2, for instance, we have 



<S\i)>= ijj(t ■ cos6,t ■ sm6) dO, 
Jo 

r2n 

< S t ,ip >— / t ■ ij)(t ■ cos9,t ■ sin0) d9, 
Jo 

and so by Proposition [3 

r-t r2n 

< B t ,ip >— / [/ u ■ ip(u • cos 9, u ■ sin 6) d6] du, 
Jo Jo 

which the reader may want to rearrange, using Fubini, into the standard 
formula for integration in polar coordinates over the disk of radius t; but 
note we have no assumptions like u t > 0". 

Note also that Bq = 0, whereas S° and B° are constants times the Dirac 
distribution at the origin (use (jHJ). The constants are the "area" of the 
unit sphere, or the "volume" of the unit ball, in the appropriate dimension. 
Explicitly, 

S° = 2-5(0), 2n-S(0), 4tt-5(0), (7) 

and 

B° = 2-5(0), n- 5(0), —-5(0) (8) 

in dimensions 1,2, and 3, respectively. 

We shall also have occasion to consider the distribution (of compact sup- 
port) t ■ S f on R 3 as well as its projection p*(t ■ S f ) on the iq/-plane (using 
functorality of T>' c with respect to the projection map p : R 3 — » R 2 ). 

We insert for reference two obvious "change of variables" equations. Re- 
call that H t : R n —> R n is the homothetic transformation "multiplying by f . 

We have, for any vector field F on R n (viewed, via principal part, as a map 
Rn Rn y 

div (F o H t ) = t ■ (div F) o H t , (9) 

and 

t n f <f>oH t = [ 0. (10) 

JBx JBt 
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4 Divergence Theorem for Unit Sphere 



The Main Theorem of vector calculus is Stokes' Theorem: Jg^uj = Ldu, 
for uj an (n — l)-form, 7 a suitable n-dimensional figure (with appropriate 
measure on it) and d'y its geometric boundary. In the synthetic context, the 
theorem holds at least for any singular cubical chain 7 : I n — > M (I n the 
n-dimensional coordinate cube), because the theorem may then be reduced 
to the fundamental theorem of calculus, which is the only way integration en- 
ters in the elementary synthetic context; measure theory not being available 
therein. For an account of Stokes' Theorem in this context, see ^T] p. 139. 
Below, we shall apply the result not only for singular cubes as in loc.cit., 
but also for singular boxes, like the usual (7 : R 2 — > R 2 ,[0,2n] x [0,1]), 
"parametrizing the unit disk B by polar coordinates" , 

7(0, r) = (rcos0,rsin0). (11) 

We shall need from vector calculus the Gauss-Ostrogradsky "Divergence The- 
orem" 

flux of F over d'j = (divergence of F), 

with F a vector field, for the geometric "figure" 7 = the unit ball in R n . For 
the case of the unit ball in R n , the reduction of the Divergence Theorem to 
Stokes' Theorem is a matter of the differential calculus of vector fields, differ- 
ential forms, inner products etc. (See e.g. [7] p. 204). For the convenience 
of the reader, we recall the case n = 2. 

Given a vector field F(x,y) = (F(x,y), G(x,y)) in R 2 , apply Stokes' 
Theorem to the differential form 

uj := — G(x, y)dx + F(x, y)dy 

for the singular rectangle 7 given by (fTTj) above. Then, using the equational 
assumptions on cos, sin and their derivatives, we have 

j*(dx) = cos Qdr — r sin 6d6 
7* (dy) = sin Qdr + r cos 6d6 
Y{dx A dy) = r (dr A d6) 

Since duj = (dG/dy + dF/dx) dx A dy = div (F) dx A dy, then 

Y(duj) = div (F) r (dr A d6) 
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On the other hand, 

7*cj = (F sm9-G cos 9) dr + (F r cos6 + Gr sin#) d9, (12) 
(all F, G, and F to be evaluated at (r cos 9, r sin 9)). Therefore 

[ du= f 1 div(F) r dr d9; 

J-y JO JO 

this is f Bi div (F) dA. On the other hand by Stokes' Theorem f du = 
JgryUi which is a curve integral of the 1-form (|12j) around the boundary of 
the rectangle [0, 27r] x [0,1]. This curve integral is a sum of four terms 
corresponding to the four sides of the rectangle. Two of these (corresponding 
to the sides 9 = and 9 = 2n) cancel, and the term corresponding to the 
side where r = vanishes because of the r in r (dr A dO), so only the side 
with r = 1 remains, and its contribution is, with the correct orientation, 

/^(F(cos 9, sin 9) cos 9 + G(cos 9, sin 9) sin 9) d9 = f F n ds 
Jo JSi 

where n is the outward unit normal of the unit circle. This expression is 
the flux of F over the unit circle, which thus equals the divergence integral 
calculated above. 



5 Time Derivatives of Expanding Spheres and 
Balls 

We now combine vector calculus with the calculus of the basic ball- and 
sphere-distributions, as introduced in Section 3, to prove the following result: 

Theorem 8 In R n (for any n ), we have, for any t, 

|S> = *.A( B <), 

(A = the Laplace operator). 

Proof. We consider the effect of the two expressions on an arbitrary function 
ip. We have 

< — S 1 *, ib >=< S,—(ibo H t ) > by various definitions 
at dt 
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=< S,u ^ {Vip(H t {u)) ■ u > by differential calculus 

= ffux over S of (Vi/j ° H ) by special property of 5* 

=< B, diviVip o if*) > by divergence Theorem 

= t<B, div(Vip) o H* > by © 

= t < B, (Aip) o H l > by definition of A 

= t < B f , A(ip) >= t < A(_B*), t/> > by various definitions 

from which the result follows. 

We collect some further information about t-derivatives of some of the 
t-parametrized distributions considered. From Proposition |7| and the Funda- 
mental Theorem of Calculus, we immediately derive 

= (13) 

In dimension 1, we have 

i(St) = A(5 t ); (14) 



dt 



for, 



wheras 



< AB t ,ij>=< B t ,ip" >= f t 4>"(t) dt, 



i-t 

and the result follows from the Fundamental Theorem of Calculus. - The 
equation |T4|) implies the following equation if n = 1; we shall prove that it 
also holds if n > 2: 

t-j t (S t ) = (n-l)S t + t-A(B t ). (15) 
For, differentiate St = t n ~ l ■ S l to get 

±(S t ) = (n-l)t"- 2 .S t + t"- l j t (S t ), 



in-2 



which by Theorem |H] and the definition of B t in terms of B l is = (n— l)t r 
S f + A(B t ). Multiplying this equation by t and using the defining equation 
St = t n_1 5* gives the result. 
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We we shall finally argue that 



' dt 




(16) 



For, differentiating the defining equation B t = t n ■ B l gives d/dtB t = nt 71 ^ 1 ■ 
B l + t n ■ djdtB 1 '. Now the left hand side here is S t , by (JT3J), so we conclude 
that t n ~ l ■ S* = nt n ~ l ■ B l + t n -d/ dtB 1 . If t were invertible, we would conclude 
by cancelling t n ~ l in this equation. But since the equation holds for all t, 
we may cancel it in any case: a consequence of the integration axiom is the 
Lavendhomme Cancellation Principle, which says that if t ■ g(t) = for all t, 
then g(t) = for all t, see jS] Ch.l Prop. 15. Applying this principle n — 1 
times then yields (fTB|) . 

6 Wave equation 

Let A denote the Laplace operator Y^9 2 /dxf on R n . We shall consider the 
wave equation (WE) in R n , (for n = 1, 2, 3), 



as a second order ordinary differential equation on the Euclidean vector space 
V c (R n ) of distributions of compact support; in other words, we are looking 
for functions 



so that for all t £ R, Q(t) = A(Q(t)) (viewing A as a map V' c (R n ) -> 



Consider a function / : R — > V, where V is a Euclidean vector space 
(we are interested in V — T>' c (R n )) . Then we call the pair of vectors in V 

consisting of /(0) and / (0) the initial state of /. We can now, for each of the 
cases n = 1, n = 3, and n = 2 describe fundamental solutions to the wave 
equations. (The case n = 2 is less explicit, and is derived "by projection" 
from the one in dimension 3.) By fundamental solutions, we mean solutions 
whose initial state is either a constant times (5(0,0)), or a constant times 



dP_ 
dt 2 



Q(t) = AQ(t) 



(17) 



Q : R -> V' c (R n ) 



(0,5(0). 
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Theorem 9 In dimension 1: The function R — > T>' C (R) given by 



t ^ S\= S t ) 



is a solution of the WE; its initial state is 2(5(0), 0). 
The function R — > T>' C (R) given by 



t i-» B, 



is a solution of the WE with initial state 2(0, 5(0)). 

Proof. We have d/dt(B t ) = S t by flI3J), and d/dt(S t ) = A(B t ), by (JUJ). This 
establishes the WE for B t . Since A and d/dt commute, it therefore follows 
that WE also holds for d/dt(B t ) = St- The initial position of the solution St 
is S = S° = 25(0), by (JTj), and the initial velocity A(B ) by (JUJ), which is 
since Bq = 0. The initial state of the solution B t is Bq = 0, and the initial 
velocity is So, as we already calculated, so is = 25(0). 

Theorem 10 In dimension 3: The function R — > TJ' C (R 3 ) given by 



is a solution of the WE with initial state 47r5(0),0). The function R — > 
V' C (R 3 ) given by 

t^t-S* 

is a solution of the WE with initial state 47r(0,5(0)). 
Proof. We calculate first d/dt of t ■ S l , using Theorem |H| 



t h+ + t 2 A(5*) 



d_ 

It 



(t-S l ) = S f + t 2 ■ A(B*) 



(18) 



and so by Theorem |H] again, 



dt 2 



(t-S l )=t- A(B l ) + 2-t- A(B*) + t 2 ■ A(— B l ) 

(JjL 

= 3-t-A(B t )+t-A(t-^B t ) 



= 3-t-A(B t )+t-A(S t -3B t ), 
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using (JSJ), and now by linearity of A, the terms involving A(B l ) cancel, so 
we are left with the equation 

^ 2 (t-S t ) = A(t.S t ) } (19) 

which establishes WE for t ■ S l . 

Since d/dt and A commute, and since t ■ S l is a solution, then so is its 
t-derivative (calculated in (fTHj) above), i.e. S* + t 2 ■ A(B t ) is a solution. The 
assertions about initial position and velocity follow from (J7J), (using Theorem 
[S] to calculate the initial velocity of the solution S t + t 2 A{B t )). 

Recall that we considered the orthogonal projection p : R 3 — > R 2 . Apply- 
ing covariant functorality, we get for any distribution Q on R 3 of compact 
support a distribution p*(Q) on R 2 , also of compact support. 

Theorem 11 In dimension 2: The function R —> T>' C (R 2 ) given by 

t^p^S* + t 2 A{B t )) 

is a fundamental solution of the WE in dimension 2; its initial state is 
4tt(<5(0), 0).The function R -> V' C (R 2 ) given by 

t^p^t-S*) 

is a fundamental solution of the WE in dimension 2; its initial state is 
4tt (0,5(0)). 

(Note: The S l and B l in the statement of the Theorem are those of R 3 .) 

Proof. The fact that the distributions in question are solutions of the 
WE is immediate from the Proposition |5] ( "p* commutes with A" ) and from 
the fact that : V C (R 3 ) — ► T>' C (R 2 ) is linear, and hence commutes with 
formation of d/dt; also, T>' c (p) sends Dirac distribution at G -R 3 to Dirac 
distribution at G R 2 , so the initial values and velocities are as claimed. 

An explicit integral expression for the two fundamental solutions here, 
obtained by projection, requires more assumptions, in particular, a square 
root formation, as is known from classical descriptions of the solutions in 
terms of "Poisson's kernel". We may express this by saying that the dis- 
tributional solutions presented exist under our weak assumptions, but that 
they are not presented by functions (densities). 
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We haven't touched the notion of support, but when defined (in a context 
where it makes sense), the two fundamental solutions , S 1 * and t-S l in dimen- 
sion 1 and 3 will have support only on the geometric sphere of radius t (which 
is of "codimension" 1), whereas the solution p*(t ■ £*) will have support in 
the direct image in R 2 of S' C i? 3 , and be of codimension 0. This accounts 
for the Huygens Principle that in a 2-dimensional world, sounds cannot be 
sharp signals, cf. e.g. [IS] p. 227. 

One might of course also derive one-dimensional fundamental solutions 
by orthogonal projection along q : R 3 — > R. Since fundamental solutions are 
unique modulo constants, we conclude that g*(S* + t 2 A(I?*)) is proportional 
to the 1-dimensional 5* (whose support is a 2-point set). 

Combining Theorem El with Proposition we can obtain information 
about S l , and other spheres and balls, for nilpotent t. As examples, we shall 
prove 

Proposition 12 Ift 5 = 0, then in dimension 1, 
and in dimension 3, 

S* = 47rM(0) + |-A(5(Q))]. 

Proof. We prove the second assertion only. (The proof of the first one is 
similar, using Theorem 0) We already observed in (J7J) that, in dimension 3, 

S° = 4vr-(5(0). 

Now the two expressions above are both solutions to WE with initial state 
(0, 47r5(0)) - the left hand side by Theorem and the right hand side by 
Proposition 0J with r = A, v = 0, w = 47r<5(0). 
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